Phonon-mediated thermal conductance of mesoscopic wires with rough edges 
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We present an analysis of acoustic-phonon propagation through long, free-standing, insulating wires 
with rough surfaces. Owing to a crossover from ballistic propagation of the lowest-frequency phonon 
mode at < oji = -nc/W to a diffusive (or even localized) behavior upon the increase of phonon 
frequency, followed by reentrance into the quasiballistic regime, the heat conductance of a wire 
acquires an intermediate tendency to saturate within the temperature range T ~ huji/ks- 



o 
Q 

> 

o 

0\ 
0\ 



X3 

o 



X 
S3 



During recent years, low-temperature heat transport 
experiments on electrical insulators have been extended 
to mesoscopic systemaH, where the wavelength of ther- 
mal phonons can be comparable to the geometrical size 
of the device. In this regime, phonon transport through a 
thermal conductor such as an electrically insulating solid 
wire, formed from an undoped semicfliiductor may ex- 
hibit ballistic waveguide propagatiorflu. This possibil- 
ity has stimulated interest in the guided-wave, phonon- 
mediated heat conductance, x{T) of ballistic wires (with 
a width W much smaller than the-|k;ngth L) connecting 
a heat reservoir to a thermal bathBB. 

In the present paper, we analyse the low-temperature 
{ksT ^ hcn/W, where c is the sound velocity) heat 
transport in relatively long {L/W ^ 100) free-standing 
insulating wires by taking into account the effect of sur- 
face roughness. The idea behind this analysis is based on 
the assumption that, in a long wire, the wire edge or sur- 
face roughness may result in strong scattering, and even 
in the localisation of acoustic waves in the intermediate- 
frequency range, whereas the low-frequency part of the 
phonon spectrum would always have ballistic properties 
due to the specifics of sound waves. In the high-frequency 
part of the spectrum, phonons would have quasiballistic 
properties, too. This may result in a non-monotonic tem- 
perature dependence of the thermal conductance of such 
a system. 

To verify the possibility of the existence of such a 
regime, in principle, we investigate the dependence on 
frequency uj of the transmission coefficient {T(lu)) using 
a simplified model of a solid waveguide which has been 
chosen to reflect two features of this problem: the influ- 
ence of roughness on the propagation of vibrations and 
the suppression of scattering on the roughness upon the 
decrease of the excitation frequency. We approach the 
problem numerically, by studying the transmission co- 
efficent averaged over many realizations of a wire char- 
acterized by a gijffin distribution of length scales in the 
surface roughnesaj and using this to find the heat conduc- 
tance, >c{T). Phonon-mediated heat transport in quasi- 
one-dimensional systems can be studied usijftg the same 
theoretical techniques as electron transportBH. However, 
in contrast to electrical conductance, which at low tem- 
peratures is determined by transport properties of elec- 
tron waves at the Fermi energy, the thermal conductance 



of a phonon waveguide is determined by all phonon en- 
ergies bj up to tiLij ~ ksT. This smears out effects of 
confinement in the transverse direction in the. temper- 
ature dependence of >c{T) in ballistic systemso, but re- 
sults in a pronounced feature in k{T) for strongly dis- 
ordered free-standing wires. The latter has the form 
of an intermediate saturation regime in >f(T) following 
a linear Tdependence at the lowest temperatures, with 
an anomalous length dependence of the saturation value, 
that scales as >Csat oc L~^l'^ for wires with a white-noise 
spectrum of roughness. 

Transmission coefficient analysis and localization of 
acoustic modes. Below, we classify phonon modes in a 
wire by the number n of nodes in the displacement am- 
plitude. The n — lowest-frequency vibrational mode 
{u! < uji = ctt/W, where c is the velocity of sound) cor- 
responds to equal displacements over the cross section 
of a free-standing wire and has a linear dispersion; oth- 
ers have frequency gaps, uJn{q) = \l {cqf + {'nnc/W)^ . 
The aforementioned feature originates from the fact that 
edge disorder suppresses the transmission of all coherent 
phonon modes at high frequencies, but has a little influ- 
ence on the n = mode at a; ^ 0, where r(a; — > 0) — > 1. 
The transmission coefficient of this mode is mainly deter- 
mined by direct backscattering, whose rate depends on 
the intensity of the surface roughness harmonic l^SW"^^ 
with wave number q uj / c. According to Rayleigh scat- 
tering theory in one dimension, the mean free path of 
this mode diverges at a; as 

Ulo) ^ (ui/iof (wy ({dW,=^/,f)) W, (1) 

even if long-wavelength Fourier components SWq are 
equally represented in the surface roughness SW{x). In 
an infinitely long wire with white-noise randomness on 
the surface, this results in a localization length for 
the lowest phonon mode which diverges at a; —> as 



(2) 



The latter statementD is based on the equivalence be- 
tweeUp. the localization problem for various types of 
wavesO. For a wire shorter than scattering yields 



1 — (r(a;)) = acj^ at cj < uJi. 



(3) 
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In contrast, at higher frequencies, all modes in a wire 
with a white-noise randomness backscatter (either via 
intra- or inter-mode process) typically at the length scale 

of I ~ W/ (^{SW/W)'^^. Hence, for frequencies u) > uji, 

the transmission coefficient tends to follow a linear fre- 
quency dependence, (r(ti;)) ~ {ujW/c) {l/L)\n{Wui/x:l, 
which is typical for diffusion in quasiballistic systemstnliJ. 
The crossover from the low-frequency regime to the 
intermediate-frequency range in a long enough wire can, 
therefore, be nonmonotonic, with a pronounced fall to- 
wards zepa at w < wi, similar to that discussed by 
Blencowell3 in relation to the phonon propagation in thin 
films. As a result, an irregular wire may exhibit ballis- 
tic phonon propagation at low frequencies, whereas, at 
higher frequencies uji < uj, surface roughness would yield 
diffusive phonon propagation, or even localization. 

The numerical simulations reported below confirm the 
above naive expectations. In these simulations, we model 
the phonons in a crystalline wire cut from a thin film 
(with the thickness much less than the wire width) as lon- 
gitudinal waves in a two-dimensional strip whose width 

/ 2\^/2 

W{x) fluctuates with rms value {{SW/Wy) = 0.1 
on a length scale longer or of order W. The effect of 
the width fluctuations consists of the scattering of acous- 
tic waves propagating along the wire. The model that 
we adopt here gives a very simplified representation of 
a real system, since we ignore the existence of torsional 
and transverse bending modes of the wire excitations, 
which are known to transfer heat in adiabatic ballistic 
constrictionaj. However, it takes into account two fea- 
tures of sound waves: their scattering and the possible 
localization by the surface roughness, and the almost bal- 
listic properties at both ultra-low and high frequencies. 
In a continuum model, these lattice vibrations are de- 
scribed by a displacement field u{x,y), which obeys 2D 
wave equation inside a wire 

w^u + c^V^M^O, (4) 

Displacements obey free boundary conditions nVw — 
at y — ±{W/2 + 6Ws{x)), where s = 1,2 indicates the 
upper and lower edge of the wire, and n stands for the 
local normal direction to the wire edge. In the simu- 
lations, we discretize Eq. (^ on a square lattice with 
about 200 sites across the wire cross section and, then, 
compute (r(w)) numerically using the transfer matrix 
method for 100 disorder realisationsE^. Our numerical 
code overcomes problems of instabilitv by QL factoris- 
ing the transfer matrices at each stepllJ. Furthermore, at 
the end of each calculation, the S-matrix is checked for 
unitarity. 

Fig. 1 shows the results of such simulations obtained 
for wires with a white-noise spectrum of roughness and 
an aspect ratio L/W = 30. Both the low- and high- 
frequency asymptotic behaviour of transmission coeffi- 
cient confirm the expected non-monotonic dependence of 
(r(a;)) for a wire of length L < L^. Calculations us- 



ing longer wires (with L ^ L^,) shows a deeper fall in 
(r(a;)) over a broader range of frequencies, since phonons 
with frequencies lo lui behave as localized vibrations, 
and their transmission coefficient becomes exponentialy 
small. The inset of Fig. 1 shows the corresponding fre- 
quency dependence of the inverse localization length of 
very long wires (with L/W ~ 1000), which is approxi- 
mately in agreement with the result of Eq. (H) . There- 
fore, the transmission of phonons through rough wires 
with white-noise roughness of the edges can be charac- 
terised by the following regimes. At low frequencies, lu < 

Lu„ where co^iL) ^ coi (w^ (^{dW.^af^^ ' VW/L, 

r ~ 1 so that phonons pass through the wire almost 
ballistically. At intermediate frequencies, cj* < o; < wi, 
phonons in the lowest mode are localized on a length scale 
shorter than the wire length, and F — > 0. Upon a further 
increase of w, modes with n ^ take part in the scatter- 
ing, the multi-mode localization length increases, and at 
uj > cL/Wl the localization length becomes longer than 
the sample length, thus r estor ing a quasi-ballistic char- 
acter to phonon transportll30. The first two regimes of 
low- frequency phonon propagation {lo < toi) through a 
wire with L ^ I can be jointly described as a function of 
a single parameter, L/Io(lu) 

{T{uj < ui)) ^p{uj/uj,); p{0) = 1, p{x > 1) - e-^. (5) 

Note that the decline of F in the vicinity oi uj ^ oji 
strongly depends on the Fourier spectrum of the rough- 
ness. To illustrate this, we analyzed the effect of rough- 
ness composed of harmonics with wave numbers q re- 
stricted to two intervals: (a) < g < n/W and 
(b) ^tt/W < q < \ti/W. The result is shown in Fig. 
2 (a) and (b) , respectively. The spectral form of the ran- 
domness is relevant, since it determines the intensity of 
Bragg-type backscattering processes. Such. processes are 
the most efficient in forming localizationtj, in which an 
incident phonon in mode n with wave number k along 
the wire axis scatters elastically to mode n' with wave 
number — fc', k' = yjk"^ + {n? — n''^){n/Wy. Therefore, 
values oi q = k + k' represented in the spectrum of SWq 
identify the regions of frequencies for which the intra- and 
inter-mode Bragg-type scattering is allowed. In Fig. 2, 
the shaded frequency intervals indicate the corresponding 
conditions for two lowest modes, n = 0, 1. 

Thermal conductance. In the regime of elastic phonon 
propagation through the wire, the heat flow Q can be 
related to the transmission coefficient of phonons through 
the wire as 

Q = ^ / -^fii^nik)vnik){fi{n,k) - f2{n,k))\tnmf , 
n.rn ^ 

where w„ = is the ID velocity of a phonon in the 
mode U!n{k), and /i(2) are equilibrium distributions of 
phonons at left (right) reservoirs. When the-temperature 
difference AT between the reservoirs is smalll3, AT <^ T, 
the thermal conductance, x — Q/AT, has the form 
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duj {hioY eyiYi{hhj /ksT) 
2tt ksT^ [ejip{hLu/kBT) ~ if 



{T{co)) (6) 



For a wire with the transmission coefficient shown in 
Fig. 2(a), where {T{u))) « 1 + cj/ivi, x{T) is plot- 
ted by the dashed-hne (1) in Fig. 3, which shows the 
crossover from hnear to quadratic temperature (depen- 
dence (at T ~ 1?! = Ghoji/kBTr^) discussed in Ref.B 



tln/eh) T + {o.7k%/h) rVi9i 



(7) 



The balhstic character of heat transport in Eq. (0) is 
reflected by the independence of >f on the sample length. 

In a wire, where the transmission coefficient sufficiently 
drops at w ~ wi, as in Fig. 1, we approximate the low- 
frequency behavior of the transmission coefficient by a 
step function, (r(ti;)) = 0{uji — uj), which yields an in- 
termediate saturation of the thermal conductance at the 
temperature T ^ -di, 
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2T/di, T<i?i; 
1, < T < diL/W. 



(8) 



even localized) behavior, with a re-entrance to the quasi- 
ballistic regime. Although the model used in this calcu- 
lation has been restricted to only one (longitudinal) exci- 
tation branch in the wire spectrum, we believe that this 
feature persists also in more realistic multi-mode models 
(which take into account torsional modes and the wire 
vibrations of other polarizations), since all lowest sound 
modes are scattered by the surface roughness with the 
rate decreasing upon the decrease of the frequency. A 
drastic difference between phonon transport properties in 
different frequency intervals results in a tendency of the 
heat conductance of a wire to saturate provisionally at 
the temperature range of T ~ hc/WkB- An intermediate 
saturation value of the wire heat conductance depends on 
the length of a wire, and, in wires with length larger than 
the scattering length of phonons with frequencies w ~ wi 
has an anomalous length dependence, >Csat oc L~^/^. 

The authors thank M.Roukes and J.Worlock for at- 
tracting our attention to this problem. This work has 
been funded in parts by EPSRC and a European Union 
TMR programme. 



The numerical result shown in Fig. 3 by a solid line is in 
a qualitative agreement with such an expectation. The 
horizontal arrow indicates the saturation value expected 
from equation (^). The upper limit in the saturation 
interval mentioned in Eq. (H) indicates the restoration 
of ballistic conditions for phonon propagation at wave- 
lengths short eaaugh to avoid wave diffraction at corru- 
gated surfacealJ'EJ. 

Theoretically, the intermediate saturation x{T) « Xsat 
at low temperatures is a more robust feature in longer 



wires of length L > Z - W^/ where even the 

lowest mode, n = is localized at frequencies a'*(L) < 

1/2 

uj<uji. Here w,(L) ~ uj^ (w^V (('5W,^o)^)) ^/W/L 
is the frequency at which the localization length of n = 
acoustic mode is comparable to the wire length. In this 
case, the saturation takes place at a lower temperature 
■i?^ ~ hw^/kBj and we find that the saturation value of the 
thermal conductance within temperature interval di > 
T > has an anomalous dependence on the sample 
length, 
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(9) 



This is an example of a more general scaling law for 
white-noise roughness; for a wire with a fractally rough 

edge, (^{SWqf'^ oc g^, one obtains Xaat oc 

In summary, our analysis of phonon propagation 
through long free-standing insulating wires with rough 
surfaces has highlighted a feature in the temperature 
dependence of the heat conductance x{T), which re- 
sults from the crossover from ballistic propagation of the 
lowest-frequency phonon mode at <C wi to diffusive (or 
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where 7 = 



of CijUij 

(? and Eij — Nij'y, and Nij is the number of nearest 
neighbours of site j from cross section i. By represent- 

ing M,. = E:"o^ 4>rxri, with = \/Wcos(^^^^) 
we prepare the equation of motion for the treatment 
by the discrete transfer-matrix method, Em ^"m^T ~ 
lEmi-^nrn'-J-Iii = ^'*r. The overlap integrals, j;^ = 
(27 - 27Cos(^))5„,„ + 7E,E±/3r'^''"^'" 



7E, XijXiiij carry information about surface rough- 
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ness through the factors (3" , which indicate the existence 
= 1) or absence {(3"' = 0) of neighbouring sites. The 
width of each wire cross-section, Wi = W + 5W\ + 5W2 
is randomly generated by introducing a random function 
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i(^X^), where Nmin and N^ax are cho- 



sen to yield the desired harmonic content, a„i arc ran- 
domly taken from the interval [—w/2,w/2], and i is the 
coordinate along the wire axis. Finally, we normalise the 
fluctuation of the width to a desired r.m.s. value, so that 
SWsii) = ^{SW^)/2 (r=(i)- <r"(i) >)/v'< r'ii^ >. 
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One can also speak of a cooling time of an overheated speci- 
men due to the phonon-mediated cooling through the wire 
into a bath, with Ttath <C To. Then, the temporal evo- 
lution of the temperature can be described by equation 
C{T)dT = -Q{T)dt, where C{T) = fn-^pYkBi^f is 
the Debye heat capacity of a specimen, and V is its vol- 
ume. The cooling time can be found by integrating this 
equation with respect to the temperature, which yeilds 
t{To) = dTC{T)IQ{T). We estimate t(To) analytically 
in two limits: i) (r(a;)) « a -|- few (with a ~ 1, 6 ~ <jJi^) 
and ii) (r((j;)) ~ ^(0;, — a;) of a weakly disordered wire and 
a strongly disordered wire with the length L S> respec- 
tively. In case i), we find approximately equal to 
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FIG. 1. Ensemble-average transmission coefficient for 100 
samples versus frequency for white-noise roughness. Inset: 
the dependence of (1 — (F)) and inverse localization length 
on the frequency for w < wi. Linear regression gave power 
dependence equal to 1.9 and 1.7, respectively. The arrows 
in the inset indicate which pair of axes corresponds to each 
curve. 



In case ii) of a long wire, we interpolate Q using 
zTiksTz) [e- - 1]-^ dz « [(T/,?.)/(l + T/i),)] !^T, 
and arrive at t(To) « 245Vpu}-^{3'&»Ti + 2T^)/efi. 
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FIG. 3. Thermal conductance versus temperature for three 
different regimes of disorder. (1) corresponds to the quasi bal- 
listic regime, (2) corresponds to the disorder regime of figure 
2(b), and (3) to the white- noise regime. The horizontal arrow 
indicates the saturation value of >t. 
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FIG. 2. Transmission coefficient versus frequency for two 
coloured noise spectra of disorder described in the text. In 
(a) the dashed line represents the transmission coefficient for 
a perfect wire. 
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